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1  Overview 


This  report  reviews  our  work  over  the  past  ~3  years  on  slow  light  nonlinear  optics.  Our  work 
consists  of  three  areas:  1)  modeling  studies  on  the  performance  limitations  of  slow  light  waveg¬ 
uides,  2)  design,  fabrication  and  characterization  of  a  novel  coupled-resonator  slow  light  waveguide 
structure  with  low  loss,  and  3)  experimental  studies  of  harmonic  generation  processes  in  plasmonic 
nanostructures. 

Of  the  most  fundamental  nature  is  our  study  of  the  limitations  on  the  performance  of  slow 
light  waveguides,  both  in  the  linear  and  nonlinear  regimes.  This  work  is  based  upon  imposing  a 
constraint  on  the  amount  of  distortion  that  a  pulse  can  experience  upon  propagation  through  a  slow 
light  waveguide.  Much  of  this  work  has  already  been  published  (and  reported  previously),  such  as 
the  effect  of  pulse  distortion  on  the  nonlinear  phase  shift,  in  that  improvement  in  nonlinear  sensi¬ 
tivity  ceases  after  a  certain  number  of  resonators  [1].  Another  study  showed  how  adjacent  pulses 
interact  in  a  nonlinear  slow  light  waveguide  [2],  where  certain  slow  light  waveguide  architectures 
(e.g.  coupled-resonator  optical  waveguides,  or  CROW)  suffer  from  relatively  strong  pulse  to  pulse 
interaction,  and  others  (e.g.  resonant-enhanced  Mach-Zehnder  interferometers,  or  REMZI)  do  not. 
We  have  also  performed  studies  in  the  linear  regime,  comparing  the  three  most  common  slow  light 
waveguide  architectures  (CROW,  REMZI,  and  side-coupled  integrated  sequence  of  spaced  optical 
resonators,  or  SCISSOR)  in  terms  of  achievable  bandwidth  delay  product  in  both  the  positive  and 
negative  (i.e.  superluminal)  group  delay  regimes  [3].  The  results  of  this  study  show  that  the  REMZI 
architecture  universally  achieves  the  greatest  bandwidth  delay  product.  This  is  due  to  the  greater 
degrees  of  freedom  available  in  the  REMZI  structure.  In  section  2,  an  overview  of  our  most  recent 
work  (to  be  completed)  on  the  ultimate  performance  of  slow  light  waveguides  will  be  discussed, 
where  the  efficiencies  of  the  three  canonical  architectures  are  compared  to  each  other  and  to  the 
ultimate  performance  predicted  by  digital  filter  optimization. 

We  have  also  been  working  on  fabricating  a  novel  photonic  crystal  based  slow  light  waveguide 
structure  in  silicon  for  use  in  nonlinear  optical  studies.  The  novel  design  reduces  the  out  of  plane 
scattering  loss  that  virtually  all  other  designs  suffer  from.  This  loss  mechanism  is  exacerbated  in 
high  index  contrast  materials  systems,  leading  to  a  tradeoff  between  bandwidth  delay  product  and 
loss.  Our  design  uses  relatively  large  cross-section  single  mode  rib  waveguides  and  deep  trench 
etching  to  simultaneously  achieve  high  index  contrast  and  very  low  loss,  as  originally  proposed 
in  [4].  We  have  developed  some  new  nanofabrication  methods  in  order  to  realize  this  structure, 
which,  along  with  characterization  results  [5],  are  discussed  in  section  3. 

Finally,  we  have  been  investigating  “slow  light”  effects  on  the  nanoscale.  In  particular,  we  have 
been  interested  in  very  small  apertures  in  metal  films,  which  serve  as  tiny  cylindrical  waveguides, 
and  other  optical  “antennae.”  When  operating  these  waveguide  just  below  the  cutoff  wavelength 
for  the  TEMqo  mode,  the  optical  density  of  states  increases,  leading  to  significantly  reduced  group 
velocity  [6]  and  enhanced  intensity  [7].  We  have  been  exploiting  this  effect  for  nonlinear  optical 
processes  such  as  second-harmonic  generation  [8]  and  third-harmonic  generation  [9].  These  results 
are  overviewed  in  section  4. 


1 


2  Performance  limitations  of  slow  light  waveguides 


Ever  since  the  demonstration  of  dramatic  reduction  in  the  group  velocity  of  light  [10],  there  has 
been  growing  interest  in  the  physics  and  applications  of  slow  light.  One  of  the  more  flexible 
implementations  of  slow  light  is  via  coupled  optical  microresontors  [11,  12,  13,  14]  (note  that 
related  structures  in  the  microwave  regime  have  been  studied  since  the  1960’s  [15,  16]).  A  particular 
application  of  interest  for  slow-light  structures  lies  in  fixed  and  tunable  optical  delay  lines  [17,  18, 
19,  20].  The  “length”  of  a  delay  line  is  often  characterized  by  the  number  of  pulse  delays  [21],  a 
metric  which  is  related  to  the  bandwidth-delay  product. 

One  of  the  important  questions  in  the  design  of  slow-light  waveguide  structures  is  the  maximum 
achievable  number  of  pulse  delays.  While  in  the  absence  of  loss,  a  long-enough  structure  could 
always  be  designed  to  meet  a  given  criterion,  a  more  interesting  question  is  the  ultimate  efficiency 
of  a  structure,  such  as  measured  by  the  number  of  pulse  delays  per  resonator.  This  efficiency  metric 
is  highly  architecture  and  design-dependent  (this  metric  also  depends  strongly  on  pulse  shape  [3]). 
We  have  taken  the  general  approach  of  discrete  time  filters  [22,  23],  where  filter  design  is  first 
performed  in  the  digital  domain  and  mapped  to  an  optical  architecture.  Discrete-time  filter  design 
concepts  [24]  have  been  widely  used  to  analyze  optical  filter  configurations,  primarily  in  terms  of 
their  linear  response.  One  of  the  limitations  of  this  approach  is  that  designs  cannot  always  be 
efficiently  mapped,  and  often,  some  constraints  are  imposed  on  the  optical  design. 

The  aim  of  this  work  is  to  compare  the  performance  of  filters  in  the  digital  and  optical  domains 
under  the  constraint  of  fixed  pulse  distortion.  The  performance  measure  for  these  filters  is  the 
number  of  pulse  delays  that  can  be  achieved,  which  is  defined  as  the  ratio  of  the  group  delay 
to  the  incident  pulse  duration.  Both  AR  (all-pole)  and  ARMA  (pole-zero)  designs  are  studied. 
The  corresponding  optical  architecture  to  the  AR  digital  filter  is  the  coupled  resonator  optical 
waveguide  or  CROW.  The  corresponding  optical  architectures  to  ARMA  digital  filters  are  side 
coupled  integrated  spaced  sequence  of  resonators  (SCISSOR)  and  ring  enhanced  Mach-Zehnder 
interferometer  (REMZI).  The  performance  (number  of  pulse  delays)  obtained  with  digital  filters 
is  better  than  the  optical  filters  because  all  the  degrees  of  freedom  for  optical  filter  design  are 
typically  not  utilized.  Often,  the  ring  length  is  a  fixed  value  in  order  to  simplify  the  mapping 
and/or  optimization  procedure,  but  that  constaint  results  in  reduced  number  of  pulse  delays.  We 
delve  into  this  issue  further  by  demonstrating  the  improvement  obtained  by  including  the  ring  length 
in  a  direct  optimization  method.  However,  due  to  to  the  increased  number  of  free  parameters  to 
optimize  with  for  a  higher  number  of  stages,  efficiency  (number  of  pulse  delays  per  resonator)  can 
drop.  The  digital  and  optical  architectures  of  AR  and  ARMA  filters  are  presented  in  section  2.1. 
In  section  2.2,  design  and  optimization  methods  for  both  digital  and  optical  filters  are  discussed. 
In  section  2.3,  the  number  of  pulse  delays  for  AR  filters  is  studied,  and  in  section  2.4,  results  for 
ARMA  designs  are  presented. 

2.1  AR  and  ARMA  filters  architectures 

The  optical  filter  architectures  and  their  transfer  matrices  are  discussed  in  this  section.  In  general, 
there  are  two  classes  of  digital  filters  -  finite-impulse  response  (FIR)  and  infinite-impulse  response 
(HR).  We  will  focus  on  IIR-type  filters,  which  are  characterized  by  poles  in  their  transfer  functions, 
and  can  be  associated  to  optical  resonaces.  The  basic  HR  filter  consists  of  all-poles  (but  may 
have  zeros  at  the  origin)  and  is  termed  autoregressive  (AR),  while  an  HR  filter  that  contains  zeros 
outside  the  origin  is  autoregressive  moving  average  (ARMA).  In  general,  ARMA  filters  have  an 
advantage  of  tailoring  the  phase  response  by  altering  the  locations  of  the  zeros  which  can  be  useful 
for  applications  like  dispersion  compensation  and  linearized  phase  modulation. 
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Figure  1  shows  the  generic  CROW  architecture  [11]  which  consists  of  a  sequence  of  N  directly- 
coupled  microring  resonators.  Each  ring  has  a  physical  length,  or  circumference,  L,  with  free 
spectral  range  FSR  =  c/nL,  where  n  =  3.0  is  the  assumed  effective  refractive  index  (Si),  and  a 
slow-light  waveguide  consists  of  N  such  coupled  resonators.  The  coupling  coefficients  are  denoted 
by  (Tj,  where  i  indexes  the  coupling  regions,  numbering  from  1  to  A/^-|-l.  We  assume  that  all  couplers 
are  ideal,  i.e.,  no  coupling  loss  and  no  backward  reflection,  so  that  operation  is  solely  determined 
by  the  field  coupling  coefficients  r  and  a  that  are  related  by  =  1 .  The  corresponding  digital 

filter  consists  of  y/N  zeros  at  the  origin  and  N  independent  poles,  which  is  termed  an  AR  lattice 
filter. 


Eout  (N  even) 


Figure  1:  CROW  architecture 


Figure  2  illustrates  a  sequence  of  waveguide-coupled  resonators,  forming  a  SCISSOR  [13,14]. 
The  unit  cell  transfer  function  for  the  SCISSOR  architecture  is  given  by  the  equation 


jj  _  Font 

J^ring  — 

-^in 


T  - 

1  —  ’ 


(1) 


where  cj)  =  kL  and  k  =  2n7r/A.  The  transfer  function  possesses  one  pole  (pi  =  1/r)  dependent 
on  its  zero  (zl  =  r),  and  is  a  constrained  ARMA  filter.  For  N  resonators,  the  transfer  function  is 


tN 


simply 


Figure  2:  SCISSOR  architecture  (top)  and  unit  cell  (bottom). 

The  architecture  of  the  REMZI  filter  with  a  single  ring  is  shown  in  Eig.  3,  where  ki  represents 
the  MZI  coupling  coefficients  and  cpR  and  4>mzi  are  the  phase  shifts  in  the  ring  and  MZI  respectively. 


3 


The  output  of  the  unit  cell  is  expressed  by  the  following  transfer  matrix  relation: 
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In  the  above  equation,  Hring  is  the  transfer  function  of  the  ring-waveguide  configuration  given  by 
equation  1.  For  N  stages,  the  corresponding  digital  filter  consists  of  N  zeros  and  N  poles  which 
are  independent  of  each  other;  this  architecture  is  represented  by  an  ARMA  lattice  filter.  The 
independent  zeros  of  the  REMZI  architecture  can  be  used  to  engineer  its  phase  response  [25,  22]. 


2.2  Filter  design  and  optimization 

This  section  presents  the  design  and  optimization  methodologies  we  have  developed  over  the  course 
of  our  research  for  digital  and  optical  filters.  Two  methods  are  studied  to  maximize  the  group  delay 
under  the  constraint  of  fixed  output  pulse  distortion.  The  first  method  is  called  the  bandwidth 
tuning  method,  where  a  prototype  linear  frequency  response  is  selected  for  the  optical  filter.  Next, 
the  optical  architecture  is  analyzed  using  the  z-transform  [26].  Then,  an  optimized  discrete  filter 
is  designed  to  give  the  same  frequency  response  as  the  prototype  response  desired  from  the  optical 
architecture.  Next,  a  mapping  algorithm  [27]  is  derived  to  synthesize  the  parameters  of  the  optical 
architecture  from  the  discrete  filter.  The  synthesized  optical  filter  is  then  simulated  using  electro¬ 
magnetic  models  and  its  linear  response  is  verified  to  be  the  same  as  that  of  the  discrete  filter. 
The  bandwidth  of  the  prototype  filter  is  tuned  in  order  to  obtain  5%  output  pulse  distortion.  The 
mapping  algorithm  between  the  digital  and  optical  domains  has  limitations  that  prevent  achieving 
the  desired  filter  characteristics  with  the  optical  architecture.  Hence,  the  group  delay  of  the  optical 
architecture  obtained  from  the  bandwidth  tuning  method  is  lower  than  that  of  the  prototype  digital 
filter. 

The  second  method,  called  the  direct  optimization  method,  involves  optimizing  the  filter  co¬ 
efficients  directly  in  the  digital  and  optical  domains.  In  the  digital  domain,  the  position  of  the 
poles  (AR,  ARMA)  and  zeros  (ARMA)  are  directly  tuned  to  maximize  the  group  delay  under  the 
constraint  of  5%  output  pulse  distortion.  In  the  optical  domain,  the  coupling  coefficients  and  phase 
shifts  of  the  optical  architecture  are  tuned  to  maximize  the  group  delay  under  the  same  distortion 
constraint.  This  method  gives  higher  group  delay  in  both  domains  for  the  same  distortion  con¬ 
straint  due  to  more  degrees  of  freedom.  In  addition,  direct  optimization  narrows  the  performance 
gap  between  the  digital  and  optical  designs. 
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Pulse  distortion  is  calculated  based  on  the  following  equation  [1] : 

„  _  f-oo 

JTjBUtWt 

where  tgd  is  the  average  group  delay  of  the  optical  architecture  and  is  calculated  based  on  the  first 
moment  of  the  output  pulse  intensity  as  given  by 

r^\EoumHt 

The  group  delay  is  maximized  for  a  10  ps  input  pulse  duration. 

2.3  Analysis  of  AR  digital  and  optical  filters 

This  section  discusses  the  design  and  evaluation  of  optical  AR  lattice  filters  using  the  digital  filter 
design  approach  and  compares  the  performance  of  AR  filters  in  the  digital  and  optical  domains. 

We  first  look  at  the  digital  AR  lattice  filter.  Figure  2.3  plots  the  number  of  pulse  delays  for 
the  two  optimization  methods.  It  should  be  noted  that  with  the  bandwidth  tuning  method,  the 
prototype  filter  response  is  maximally  flat  with  unit  magnitude  over  the  design  bandwidth.  The 
design  from  the  bandwidth  tuning  method  is  used  as  the  starting  point  for  direct  optimization.  It 
is  clear  that  direct  optimization  produces  a  significant  increase  in  efficiency,  approaching  0.5  pulse 
delays/stage. 


Figure  4:  Number  of  pulse  delays  in  the  digital  domain  (AR)  versus  number  of  stages  constrained 
by  5%  output  pulse  distortion.  Digital  designs  from  direct  optimization  and  bandwidth  tuning 
methods. 

In  order  to  implement  the  bandwidth  tuning  method  for  CROW,  we  use  a  procedure  developed 
by  Melloni  [28],  where  a  prototype  filter  is  synthesized  using  a  classic  microwave  filter  design  and 
then  substituted  with  optical  components,  such  as  coupled  ring  resonators.  Direct  optimization  is 
then  based  upon  using  this  design  as  an  initial  guess.  The  second  design  method  shows  an  increase 
in  the  number  of  pulse  delays  per  stage,  as  shown  in  Figure  5.  It  should  be  noted  that  the  mapping 
algorithm  is  not  efficient  for  arbitrary  pole  locations  and  when  the  poles  get  close  to  the  unit  circle, 
thus,  failing  to  preserve  the  digital  filter  characteristics  in  the  optical  domain. 
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Figure  5:  Number  of  pulse  delays  (AR)  versus  number  of  stages  for  optical  designs  with  5%  output 
pulse  distortion  from  direct  optimization  and  bandwidth  tuning  methods. 


In  order  to  compare  the  direct  optimization  designs  in  the  digital  and  optical  domains,  the 
number  of  pulse  delays  is  plotted  versus  the  number  of  stages  in  Figure  6.  The  discrepancy  in  the 
number  of  pulse  delays  between  the  two  domains  increases  for  higher  number  of  stages.  This  is 
mainly  due  to  the  fact  that  all  of  the  degrees  of  freedom  are  not  used  in  the  optimization  for  the 
CROW  designs.  The  additional  degree  of  freedom  is  the  ring  length.  By  including  variable  ring 
length  in  the  direct  optimizaion,  the  performance  of  the  optical  design  nearly  matches  that  of  the 
digital  design. 


Figure  6:  Number  of  pulse  delays  (AR)  versus  number  of  stages  with  5%  output  pulse  distortion. 
Digital  and  optical  designs  from  the  direct  optimization  method. 


2.4  Analysis  of  ARMA  digital  and  optical  filters 

ARMA  architectures  have  poles  and  zeros,  increasesing  their  design  flexibility.  In  general,  the 
performance  of  ARMA  filters  in  the  digital  and  optical  domains  is  better  than  that  of  AR  designs. 


6 


We  again  use  the  two  design  methods  for  the  ARMA  filters.  In  the  first  method,  we  tune  the  filter 
bandwidth  which  moves  the  pole-zero  locations  in  the  digital  domain,  and  in  turn  determines  the 
coupling  coefficients  in  the  optical  domain.  The  number  of  pulse  delays  is  plotted  in  Figure  2.4  for 
the  bandwidth  tuning  method.  For  N  =  1,  SCISSOR  and  REMZI  designs  have  the  same  group 
delay  for  5%  pulse  distortion  as  a  single  stage  of  REMZI  reduces  to  a  single  waveguide  coupled 
resonator. 


Eigure  7;  Number  of  pulse  delays  (ARMA)  versus  number  of  stages.  5%  output  pulse  distortion. 
Digital  and  optical  designs  from  direct  optimization  method. 

With  the  direct  optimization  method,  the  design  from  the  bandwidth  tuning  method  is  used 
as  a  starting  point  for  subsequent  optimization.  In  the  digital  domain,  direct  optimization  of 
the  locations  of  the  poles  and  zeros  is  performed  to  maximize  the  group  delay  under  the  same 
5%  distortion  constraint.  Eor  the  optical  architectures,  the  coupling  coefficients  (SCISSOR  and 
REMZI)  and  phase  shifts  (REMZI)  are  directly  optimized.  Similar  to  the  AR  designs,  the  direct 
optimization  method  results  in  improved  ARMA  designs.  In  Figure  2.4,  we  plot  the  number  of 
pulse  delays  for  both  digital  and  optical  designs  obtained  from  the  direct  optimization  method. 
The  results  for  REMZI  are  plotted  for  both  fixed  and  variable  ring  lengths.  Variable  ring  length 
results  are  not  plotted  for  SCISSOR  as  it  did  not  produce  any  significant  improvement.  It  can  be 
observed  that,  by  including  ring  length  in  optimization,  the  performance  improves  and  gets  closer 
to  digital  domain.  However,  there  is  still  discrepancy  due  to  the  difficulty  in  optimization  for  higher 
number  of  stages  due  to  increased  free  parameters  to  optimize  with. 

The  transfer  function  for  digital  design  with  V  =  3  is  presented  in  Figure  8  for  5%  output  pulse 
distortion  level  under  direct  optimization.  The  transfer  functions  for  the  optical  ARMA  designs  for 
N  =  3  are  presented  in  Figure  9  (REMZI  and  SCISSOR)  along  with  the  corresponding  pole  zero 
diagrams.  The  output  pulse  profiles  for  the  ARMA  designs  is  presented  in  Eig.  11 

2.5  Conclusion 

Digital  filter  design  concepts  were  applied  to  study  the  limits  on  the  performance  of  slow-light  optical 
filters  in  terms  of  the  number  of  attainble  pulse  delays  per  stage,  as  well  as  in  numerous  other  studies 
based  upon  nonlinear  pulse  propagation.  Traditional  design  methods  based  upon  prototype  digital 
filters  do  not  produce  optimal  performance  in  either  domain,  and  there  are  discrepencies  between 
the  two  domains  due  mainly  to  limitations/constraints  in  the  procedure  that  maps  poles  and  zeros 
into  optical  coupling  and  phase  shift  coefficients.  However,  direct  optimization  improved  the  filter 
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Figure  8:  Transfer  function  for  the  digital  ARM  A  filter  under  direct  optimization  -  N  =  3  and  5% 
output  pulse  distortion. 


Figure  9;  Transfer  functions  (top)  and  pole  zero  diagrams  (bottom)  for  the  REMZI  and  SCISSOR 
optical  architectures  under  direct  optization  -  N  =  1  and  5%  output  pulse  distortion. 


Figure  10:  Output  pulse  profiles  for  the  ARMA  designs  under  direct  optimization.  Digital  (top) 
and  optical  (SCISSOR  and  REMZI,  bottom)  with  A"  =  1  and  5%  pulse  distortion. 
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performance  in  both  domains.  We  showed  that  the  primary  cause  of  the  discrepency  between 
digital  and  optical  designs  is  due  to  constraining  the  optical  designs  to  a  fixed  ring  length.  When 
the  ring  length  was  also  included  as  a  free  parameter  in  optimization,  the  performance  improved 
and  became  closer  to  the  corresponding  digital  filter  designs.  From  our  results,  we  estimate  that 
the  CROW  optical  architecture  can  produce  0.5  pulse  delays/stage,  while  SCISSOR  and  REMZI 
produce  0.7  and  1.5,  respectively. 

Our  work  has  consistently  demonstrated  that  the  REMZI  slow-light  architecture  is  the  most 
efficient  across  all  metrics  related  to  linear  and  nonlinear  effects.  However,  REMZI  architectures 
are  the  most  difficult  to  fabricate. 
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3  Low  loss  photonic  crystal  slow  light  waveguides 


We  have  completed  the  fabrication  and  testing  of  the  photonic  crystal  rib  waveguide  design.  The 
nature  of  the  design  significantly  reduces  out  of  plane  scattering  (which  occurs  via  diffraction  from 
unguided  regions  of  the  photonic  crystal)  by  using  a  large  cross-section  rib  waveguide  for  transverse 
confinement,  which  supports  single  mode  operation  with  large  mode  size.  The  major  disadvantage 
of  this  approach  is  that,  in  order  to  etch  air  gaps,  or  trenches,  into  the  rib  waveguide  to  define  the 
photonic  crystal  in  the  longitudinal  direction,  high-aspect  ratio  etching  is  required.  As  shown  in 
Figure  11,  the  thickness  of  the  Si  layer  is  5.3  /rm.  For  a  quarter-wavelegnth  air  spaced  mirror  period, 
the  width  of  the  etched  trench  must  be  ~330  nm  at  1550  nm  operating  wavelength,  requiring  an 
aspect  ratio  of  16:1  or  greater. 


Figure  11:  Micrographs  of  the  slow  light  waveguide,  a)  Cross-section  of  a  rib  waveguide  showing 
the  top  and  bottom  Si02  claddings.  The  rib  waveguide  has  a  height  oi  H  =  5.3  /rm,  a  rib  width 
of  w  =  3.8  //m,  and  a  slab  height  h.  b)  A  magnified  view  of  the  surface  of  the  high  aspect  ratio 
trenches.  Three  trenches  are  use  to  separate  each  resonator,  c)  The  trenches  shown  are  etched 
through  the  rib  waveguide,  d)  A  structure  with  25  coupled  resonators.  The  input  and  output  rib 
waveguides  can  be  seen  leading  to  the  CROW. 

The  slow-light  waveguides  are  fabricated  on  a  100  mm  SOI  wafer  using  a  series  of  surface 
micromachining  processes.  The  flow  diagram  for  the  fabrication  is  shown  in  Figure  12.  The  first 
fabrication  step  is  to  perform  optical  lithography  to  create  alignment  marks  for  the  trenches  and 
waveguides.  The  alignment  marks  are  etched  into  the  wafer  using  reactive  ion  etching  (RIE)  with 
SFe  and  O2  gases.  Next,  a  thin  layer  of  Cr  is  sputtered  over  the  wafer  followed  by  a  spin-on  layer 
of  resist  (ZEP  520A,  Zeon  Corp.).  The  resist  is  patterned  using  e-beam  lithography  (NovaNano, 
FEI)  and  then  developed.  The  pattern  is  transferred  to  the  Cr  using  RIE  with  CI2  gas.  The  Si  is 
then  patterned  using  cryogenic  deep  reactive  ion  etching  (DRIE)  using  SEg  and  O2  gases.  The  etch 
recipe  gives  smooth  vertical  sidewalls.  The  residual  Cr  mask  is  removed  using  wet  etching  (Cr-14 
etchant). 

Next,  the  rib  waveguide  is  formed  over  the  trenches  using  optical  lithography  and  another 
DRIE  step.  This  etching  step  uses  a  much  higher  pressure  than  the  previous  DRIE  step  to  create 
a  short  mean  free  path  between  reactive  gas  molecules  which  prevents  etching  at  the  bottom  of  the 
trenches.  The  waveguide  must  be  aligned  perpendicularly  to  the  trenches  to  avoid  deflecting  light 
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Figure  12:  Process  steps  used  to  fabricate  the  waveguides,  a)  An  SOI  wafer  is  processed  starting 
with  b)  an  optical  lithography  step  using  a  positive  photoresist  (Shipley  1813)  and  c)  dry  etching 
using  RIE.  d)  A  layer  of  Cr  is  sputtered  on  the  wafer  surface  followed  by  e)  a  spin-on  layer  of 
e-beam  resist.  The  resist  is  then  patterned  and  f)  the  Cr  is  etched  using  RIE.  The  Cr  mask  is  then 
used  as  a  mask  during  g)  etching  of  the  trenches,  h)  Another  optical  lithography  step  is  used  to 
create  an  etch  mask  for  the  rib  waveguide,  which  is  i)  etched  using  a  short  DRIE  step. 


out  of  the  waveguide.  The  alignment  is  better  than  a  half  degree  from  90°,  similar  to  state-of-the-art 
fabrication  quality  achieved  by  others. 

The  trenches  are  back  filled  by  wet  oxidation.  The  wet  oxidation  is  self-arresting  as  the  oxide 
from  both  sides  of  the  trenches  meet.  The  purpose  of  the  oxidation  is  to  reduce  surface  roughness 
inside  the  trenches  which  lowers  the  waveguide  loss.  Roughness  scattering  is  particularly  detri¬ 
mental  to  slow  light  waveguides  because  the  scattering  loss  scales  with  the  square  of  the  electric 
field. 

We’ve  achieved  high  aspect  ratio  etching  using  an  inductively  coupled  plasma  (ICP)  reactive 
ion  etching  (RIE)  system.  Examples  of  this  high  aspect  ratio  etching  are  shown  in  Eigure  13.  Two 
types  of  structures  have  been  made.  The  first  is  with  the  photonic  crystal  region  defined  by  air  gaps. 
In  order  to  obtain  the  advantages  of  reduced  out  of  plane  scattering,  vertical  symmetry  should  be 
maintained,  which  is  achieved  by  using  chemical  vapor  deposition  (PECVD)  of  an  oxide  cladding 
layer.  While  there  is  some  deposition  into  the  air  trenches  (seen  by  tapering  of  the  top  portions  of 
the  trenches),  the  air  trenches  remain  intact.  The  other  structure  has  trenches  filled  in  with  oxide, 
which  is  accomplished  by  wet  thermal  oxidation  of  the  structure  after  the  air  trenches  have  been 
etched.  This  requires  compensation  of  the  initial  trench  width  due  to  the  take-up  of  Si  material 
during  oxidation;  PECVD  oxide  is  then  deposited  as  a  top  cladding  layer  to  maintain  symmetry. 

Transmission  measurements  through  these  PCRWs  have  been  made  using  a  tunable  telecom 
laser.  The  transmission  spectra  are  then  normalized  to  the  transmission  measurement  through 
identical  rib  waveguides  (but  without  the  PC  structure).  Some  of  the  results  are  shown  in  Eig¬ 
ure  14  for  structures  with  5,  10,  and  25  resonators.  A  few  features  can  be  noticed  from  these 
results.  In  general,  the  entire  bandgap  of  the  structures  is  wider  than  the  tuning  range  of  our  laser 
(approximately  180  nm),  where  only  the  high  frequency  bandedge  and  central  defect  transmission 
resonance  are  within  range;  the  entire  bandgap  region  is  expected  to  cover  approximately  200  nm, 
corresponding  to  13%  fractional  bandwidth  (for  an  air  trench  structure,  the  width  will  be  even 
greater).  Another  feature  is  the  large  dynamic  range  between  the  transmission  windows  and  the 
bandgap  region,  exceeding  10^  for  25  resonators.  The  number  of  peaks  in  the  miniband  corresponds 
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Figure  13:  Cross-section  images  of  the  photonic  crystal  region  of  PCRWs  using  an  Si02  cover  layer 
with  air  trenches  (left)  or  using  wet  oxidation  to  fill-in  the  trenches  (right). 


to  the  number  of  resonators;  as  the  number  of  resonators  increases,  the  transmission  band  becomes 
flatter  with  sharper  roll-off,  as  expected  for  N  ^  oo.  Finally,  what  is  most  surprising  is  the  high 
transmission  through  the  PCRW,  where  up  to  70%  transmission  is  obtained  through  25  resonators. 

We  have  calculated  the  slowing  ratio  and  propagation  losses  for  these  structures.  The  slowing 
ratio  is  29  and  the  loss  is  just  under  1.5  dB/cm,  where  the  effective  path  length  has  been  used.  For 
slow  light  waveguides  to  be  used  in  all-optical  networks  as  time-delay  buffers,  they  must  be  capable 
of  significantly  reducing  the  group  velocity.  In  addition,  they  must  also  have  a  wide  bandwidth  to 
accommodate  optical  signals  and  avoid  pulse  distortion.  Therefore,  the  delay-bandwidth  product 
(DBP)  is  an  important  figure  of  merit.  For  the  25  resonator  waveguide,  the  bandwidth  over  which 
ng  varies  by  ±10%  about  the  average  Ug  is  16  nm.  The  calculated  delay  time  is  5.4  ps,  giving 
a  DBP  of  10.8.  This  is  a  relatively  large  value  and  compares  well  to  state-of-the-art  slow  light 
waveguides  (see  Table  1). 

The  DBP  could  be  increased  by  coupling  more  resonators  together.  However,  as  more  cavities 
are  added,  loss  also  increases,  limiting  the  practical  length.  Therefore  another  important  figure 
of  merit  is  the  DBP  per  unit  loss.  For  our  structure  with  25  resonators  the  DBP/loss  is  38  per 
dB,  which  is  also  shown  in  Table  1  along  with  the  values  for  several  other  slow  light  waveguides 
from  the  literature.  On  a  DBP/loss  basis,  the  waveguide  presented  here  essentially  ties  for  best 
performance.  The  normalized  DBP  (nDBP)  is  also  an  important  figure  of  merit  for  comparing  the 
performance  of  various  waveguides.  The  nDBP  is  defined  at  ng{Au}/u}),  where  uj  is  the  frequency 
and  Au!  the  bandwidth.  For  the  structure  presented  here,  the  value  is  0.31,  which  is  higher  than 
any  other  slow  light  structure. 

Despite  the  relatively  large  height  of  the  rib  waveguide,  the  overall  footprint  of  the  structure  is 
still  quite  small.  The  25  resonator  CROW  as  fabricated  occupies  a  surface  area  of  approximately 
3  X  10“^  mm^,  which  is  nearly  half  the  footprint  of  25  coupled  ring  resonators,  assuming  8  fj,m  radii 
and  200  nm  gaps.  The  rib  waveguide  CROW  footprint  could  be  reduced  further  by  shortening  the 
length  of  the  trenches.  Our  BPM  simulations  show  that  the  trenches  only  need  to  be  6  times  the 
waveguides  rib  width  w.  Therefore,  the  footprint  of  the  25  resonator  structure  could  be  reduced 
to  1.4  X  10^  mm^.  This  compact  footprint  makes  the  rib  waveguide  CROW  promising  for  on-chip 
integration  with  other  optical  structures. 
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Figure  14:  Normalized  transmission  measurements  through  PCRWs  of  5  (top),  10  (middle)  and  25 
(bottom)  defect  resonators. 


Table  I.  Comparison  of  slow  light  waveguides  from  the  literature. 


Structure 

Measurement 

(iroup 

DBF 

Propagation 

DBP/loss 

Reference 

IXE? 

Index 

Loss  (dB/cm) 

(dB-‘) 

1 .  Coupled  Rib 

Waveguide 

Resonators" 

Specual 

30.1 

10.8 

52.7 

38.0 

Presented 

Here 

2.  Coupled  PC 
Resonators 

Temporal  and 
Socctral 

129 

1 1  kf 

14.5‘ 

39.9 

Notomi  [10) 

3.  Coupled  PC 
Resonators 

Spectral 

105 

'ia‘ 

3x10' 

0.4 

Jdgerska  [28] 

4.  Coupled  Ring 
Resonators 

Temporal  and 
Spectral 

36.3" 

55" 

126" 

2.4 

Xia  [9l 

5.  PC  Band  Edge 

Spectral 

34 

78' 

150 

10.4 

Frandsen  [30j 

6.  PC  Band  Edge 

Spectral 

25 

30 

244 

7.7 

Settle  |3 1] 

7.  PC  Band  Edge 

Spectral 

42 

57 

- 

- 

Baba  [32] 

8.  PC  Band  Edge 

Spectral 

44 

15.8" 

- 

- 

ui29| 

“  Measured  for  a  25  resonator  structure  with  the  spectrum  shown  in  Fig.  6. 

*  Calculated  using  bandwidth  of  1 .2  nm.  N  =  6()  resonators.  V|  =  7.7x1 0’^c,  and  75  ps  delay  time. 

'  Estimated  loss. 

Calculated  from  paper. 

*  Calculated  for  1(X)  coupled  rings  in  CROW  configuration,  i.e.  9  pm  radii  and  2(X)  nm  gap.  The  linear  length 
is  used  to  calculate  the  propagation  loss  and  group  index,  not  the  total  path  length.  A  2  nm  bandw'idih 
centered  at  1 550  nm  is  used  to  calculate  DBF. 

^Calculated  with  1 1  nm  bandwidlli.  i.e.  bandwidth  over  which  ng  is  reported  to  be  constant. 


Figure  15:  Comparison  of  slow  light  waveguides  from  the  literature.  The  citations  can  be  found  in 
the  Optics  Express  2010  paper. 
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4  Third  harmonic  generation  from  sub- wavelength  apertures 


Extraordinary  optical  transmission  (EOT)  through  metal  films  modulated  with  a  2-D  array  of 
sub-wavelength  apertures  [29]  was  reported  in  1998.  It  is  now  well  established  that  upon  the 
condition  of  EOT,  intensity  buildup  within  the  apertures  can  occur  [30,  31],  motivating  the  study 
of  nonlinear  processes.  Indeed,  resonant  enhancement  of  second-harmonic  generation  (SHG)  by  a 
factor  of  10^  has  been  demonstrated  for  a  single  sub-wavelength  aperture  surrounded  by  periodic 
annular  corrugation  [32].  SHG  has  also  been  studied  in  arrays  of  sub- wavelength  apertures  of 
various  shapes,  using  disordered  [33,  34,  8]  and  periodic  [34,  35,  36,  37,  38]  arrangements.  The 
effect  of  the  symmetry  of  the  aperture  has  been  shown  [34],  in  that,  at  normal  incidence,  apertures 
with  inversion  symmetry  produce  much  weaker  SH  than  non-centrosymmetric  apertures,  while  at 
off-normal  incidence,  SH  can  be  produced  with  centrosymmetric  apertures  [34,  38,  8]. 

So  far,  nonlinear  optical  studies  have  been  directed  towards  effects,  although  multi-photon 
luminescence  has  been  observed  [37].  Glosely  related  are  nonlinear  optical  studies  using  bow-tie 
antennae,  where  strong  field  enhancement  is  obtained  within  the  feedgap  [39].  Recently,  high- 
harmonic  generation  was  demonstrated  utilizing  the  strong  field  enhancement  within  the  gaps  of 
bow-tie  antennae  [40].  Our  focus  is  on  third-harmonic  generation  (THG)  in  transmission  from  sub¬ 
wavelength  metallic  apertures  arrayed  in  gold  films,  where  the  effects  of  aperture  size  and  spacing 
on  conversion  are  demonstrated. 

4.1  Experimental  methods 

Arrangements  of  sub- wavelength  apertures  were  produced  in  100  nm  thick  gold  films  using  electron 
beam  lithography  (EBL).  Briefly,  a  5  nm  chromium  or  Ti02  adhesion  layer  was  sputter  deposited 
onto  the  quartz  substrate,  followed  by  100  nm  of  gold  and  20  nm  of  chromium.  ZEP520A  e-beam 
resist  of  about  300  nm  thickness  was  spin  coated.  Following  e-beam  exposure,  the  upper  chromium 
layer  was  dry  etched  with  chlorine,  and  the  e-beam  resist  removed.  The  chromium  layer  served 
as  a  hard  mask  for  argon  ion  milling  of  the  gold.  A  wet  etch  removed  the  upper  chromium  layer 
(and  likely  resulted  in  some  undercut  in  the  underlying  chromium  adhesion  layer).  For  samples 
with  Ti02  as  the  adhesion  layer,  we  used  reactive  sputtering  to  produce  the  ~5  nm  Ti02  layer  by 
introducing  O2  during  the  sputtering  process;  pure  O2  was  flowed  through  the  chamber  to  fully 
oxidize  the  thin  film  before  the  chamber  was  pumped  down  to  sputter  the  gold  film. 

The  experimental  setup  is  similar  to  the  one  used  previously  [8],  as  shown  in  Fig.  16.  A 
Ti:Sapphire  laser  is  used  at  800  nm  wavelength  and  roughly  30  fs  pulse  duration.  The  setup  allows 
the  rotation  of  the  sample  with  respect  to  the  incident  fundamental  beam  as  well  as  rotation  of 
the  detector  around  the  sample  so  that  the  radiation  pattern  can  be  measured  in  transmission. 
The  detector  is  a  blue-sensitized  PMT  (H5784-03)  and  two  spectral  filters  are  used  to  block  the 
transmitted  fundamental  at  800  nm  and  minimize  the  influence  of  broadband  background  lumines¬ 
cence  [37].  The  angular  acceptance  of  light  collection  is  roughly  1  degree  through  a  slit  of  about 
2  mm  width.  Lock-in  detection  is  performed  by  modulating  the  86  MHz  pulse  train  at  2  kHz,  and 
average  incident  power  is  75  mW  (measured  after  the  chopper),  except  where  noted. 

4.2  Results 

One  significant  difference  between  SHG  and  THG  is  that  THG  does  not  have  the  symmetry  limi¬ 
tations  of  SHG.  This  can  be  seen  in  Fig.  17  which  plots  the  fundamental  transmission,  SH  output, 
and  TH  output  versus  incidence  angle  for  a  sample  with  Gr  adhesion  layer,  and  885  nm  spacing 
of  round  holes.  Because  of  centro-symmetry  of  the  sample  at  normal  incidence,  there  can  be  no 
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Figure  16:  Experimental  setup  for  third  harmonic  generation  measurements.  The  output  of  the 
TiiSapphire  laser  passes  through  a  neutral  density  hlter  and  a  spectral  longpass  filter  (690  nm) 
before  being  focused  onto  the  sample  with  a  20  ^m.  spot  size  by  a  10  cm  focal  length  lens;  all 
patterned  areas  were  larger  than  the  spot  size  (EBL  samples  are  80  /rm  by  80  /rm).  Emission  from 
the  sample  passes  through  a  collection  lens  of  5  cm  focal  length  made  from  fused  silica,  a  spectral 
UG5  glass  hlter  (Thorlabs)  with  passband  around  225  nm  to  400  nm  to  suppress  the  fundamental, 
an  Hg  line  interference  hlter  (CVI,  265  nm  center,  25  nm  passband)  to  isolate  the  TH,  and  a  2  mm 
slit,  and  is  detected  with  the  PMT. 


SH  output  [34,  8].  However,  the  response  always  exists.  As  shown  in  Fig.  17,  the  TH  output 
peaks  with  maxima  in  fundamental  transmission,  indicating  that  the  greatest  intensity  enhance¬ 
ment  within  the  aperture  occurs  at  the  transmission  maxima,  as  expected. 

Power  scaling  of  the  SH  and  TH  signals  are  also  shown  in  Fig.  17.  The  measured  data  points 
are  plotted  on  a  log-log  scale  and  ht  to  a  linear  equation  with  slopes  of  2.0  and  2.8,  respectively, 
for  SH  and  TH,  which  represent  the  power  scaling  factors.  From  the  measurements,  it  is  clear  that 
the  samples  can  withstand  average  incident  power  levels  where  the  TH  signal  is  about  the  same  as 
the  maximum  SH  (under  strong  symmetry- breaking  conditions),  not  correcting  for  differences  in 
collection  efficiencies  and  PMT  responsivity.  Whether  this  holds  for  even  higher  harmonics  is  not 
known  yet. 

Double-angle  scans  for  SHG  and  THG  using  the  sample  with  Cr  adhesion  layer  are  shown  in 
Fig.  18,  where  both  the  incidence  and  detection  angles  are  varied.  Again,  the  difference  between  the 
SHG  and  THG  mechanisms  is  observed  in  that  there  is  minimal  SH  under  conditions  of  inversion 
symmetry  (normal  incidence  and  normal  detection),  other  than  a  small  two-photon  luminescence 
background.  The  peaks  in  the  detection  scan  correspond  to  coherent  emissions  from  the  apertures 
that  satisfy  the  following  momentum  matching  condition  [8]: 

=  nk^  +  mK,  (5) 

where  n  is  the  harmonic  order,  K  is  a  reciprocal  lattice  vector  with  |K|  =  2tiI K,  A  is  the  aperture 
spacing,  m  is  the  diffraction  order,  and  k*  represents  a  transverse  light  wave-vector.  For  a  square 
lattice,  and  assuming  that  the  optical  wavevectors  have  only  the  x  transverse  component  (since  the 
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Incident  Power  (mW) 


Figure  17;  Fundamental  transmission,  second-harmonic  output,  and  third-harmonic  output  versus 
incidence  angle  for  a  sample  of  200  nm  holes  with  885  nm  spacing  in  a  square  lattice  (left).  The 
adhesion  layer  is  5  nm  of  Cr.  Power  scaling  of  second-harmonic  and  third-harmonic  outputs  (right). 
For  THG  measurement,  the  sample  was  set  for  a  12°  incidence  angle,  while  for  SHG  measurement, 
the  sample  was  set  for  46°  incidence. 


scans  are  performed  along  only  one  axis),  this  equation  can  be  written 

siny  =  sin0 -|- (6) 
nA 

where  7  is  the  detection  angle  and  9  is  the  incidence  angle.  Therefore,  the  angular  spacing  between 
peaks  in  the  TH  measurements  are  narrower  due  to  the  shorter  wavelength  of  emission  (n  =  3).  At 
detection  angles  that  lie  in-between  these  peaks,  luminescence  background  can  be  observed.  Since 
two-  and  three-photon  luminescence  are  incoherent  emissions,  they  produce  broad  angular  spectra, 
which  is  clearly  evident.  Evidence  that  the  luminescence  emits  from  the  apertures  is  given  by  the 
fact  that  the  luminescence  is  strongest  at  the  same  incidence  angles  that  the  SH  and  TH  peaks 
occur. 


Figure  18:  SH  (left)  and  TH  (right)  emission  for  a  sample  with  Gr  adhesion  layer,  885  nm  pitch, 
and  200  nm  aperture  diameter.  The  SH  and  TH  signal  intensities  are  normalized  to  the  same 
value.  Trajectories  of  the  diffraction  peaks  follow  equation  6,  as  indicated  by  the  dashed  lines  on 
the  contour  plots. 

A  practical  issue  not  always  addressed  for  plasmonic  structures  is  the  adhesion  layer.  Often- 
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times,  a  thin  layer  of  Cr  is  used  to  promote  the  adhesion  of  Au  to  a  glass  substrate.  However, 
Cr  is  a  highly  lossy  material,  and  can  cause  significant  attenuation  of  SPP  propagation  at  the 
metal/substrate  interface  [41,  42].  We  performed  all  subsequent  TH  measurements  on  samples 
with  Ti02  adhesion  layers.  The  effects  of  different  adhesion  layers  will  be  discussed  in  a  separate 
paper. 

4.2.1  Effects  of  aperture  spacing 

Keeping  the  aperture  size  hxed  at  250  nm,  we  investigated  the  effects  of  varying  the  spacing.  The 
results  are  shown  in  Figure  19  for  both  the  transmission  of  the  fundamental  beam  and  emission  at 
the  TH  wavelength.  In  these  measurements,  the  sample  was  rotated  with  respect  to  the  incident 
beam  (varying  9)  and  detection  was  performed  at  the  zeroth  order  transmission  (7  =  —9). 
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Figure  19:  Fundamental  transmission  (left)  and  third-harmonic  output  (right)  for  samples  of  250  nm 
apertures  with  variable  spacing  in  a  square  lattice.  The  adhesion  layer  is  5  nm  of  Ti02. 


Two  effects  of  changing  the  aperture  spacing  are  readily  observable.  First,  there  is  a  shift 
to  higher  incidence  angles  at  which  the  transmission/TH  emission  peak.  Second,  there  is  a  clear 
increase  in  TH  output  for  spacings  of  885  nm  and  longer;  note  that  this  is  repeatable  across  multiple 
measurements  of  these  patterns.  Even  though  the  fundamental  transmission  peaks  also  increase 
with  aperture  spacing,  the  effect  is  not  as  dramatic. 


4.2.2  Effects  of  aperture  size 

Keeping  the  aperture  spacing  fixed  at  about  885  nm,  we  investigated  the  effects  of  varying  the  aper¬ 
ture  size.  The  results  are  shown  in  Figure  20.  Again,  there  are  two  effects  which  are  immediately 
noticed.  The  first  is  that  the  transmission  of  the  fundamental  increases  with  increasing  aperture 
size,  which  is  expected;  however,  the  second  effect  is  unexpected.  At  the  first  two  peaks  in  TH 
output  (near  zero  and  10°),  TH  decreases  with  aperture  size,  while  at  the  third  peak  (near  30°), 
TH  increases  with  aperture  size;  at  the  fourth  peak,  it  decreases  again.  Individual  apertures  have 
cutoff  frequencies,  below  which  group  velocity  is  minimized,  resulting  in  increase  in  intra-aperture 
intensity  [6,  37].  This  suggests  an  explanation  for  the  behavior  seen  at  0°,  10°,  and  50°,  where  the 
decreasing  aperture  size  should  result  in  increasing  TH.  However,  this  doesn’t  explain  the  results 
at  30°.  The  fundamental  difference  between  these  two  cases  lies  in  the  localization  of  the  enhanced 
electromagnetic  field. 

In  order  to  investigate  the  origin  of  the  scaling  of  TH  with  aperture  size,  we  performed  calcu¬ 
lations  using  Lumerical  FDTD  Solutions.  Figure  21  shows  the  results  of  these  calculations,  where 
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Figure  20:  Fundamental  transmission  (left)  and  third-harmonic  output  (right)  for  samples  of  vari¬ 
able  hole  size  with  885  nm  spacing  in  a  square  lattice.  The  adhesion  layer  is  5  nm  of  Ti02. 


transmission  of  the  funamental  is  plotted  versus  incidence  angle  for  different  aperture  sizes.  By 
comparing  with  Fig.  20,  good  agreement  with  experiment  is  obtained.  The  origin  of  the  emitted 
TH  signal  is  from  the  interior  metal  surface  of  the  aperture  and  the  substrate  at  the  bottom  of  the 
aperture,  where  E^{uj)  at  these  locations  is  the  source  polarization.  Figure  21  also  plots  the  average 
value  of  within  a  volume  that  lies  10  nm  from  the  interior  surface  of  the  aperture.  As  shown, 
the  intensity  increases  with  aperture  size  at  25°,  but  decreases  at  10°,  supporting  the  experimental 
measurements  well. 
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Figure  21:  Calculated  fundamental  transmission  versus  incidence  angle  for  different  aperture  di¬ 
ameters  (left)  with  885  nm  spacing.  Calculated  average  \E\^  within  a  10  nm  thick  bucket-shaped 
volume  versus  aperture  size  for  angles  of  incidence  corresponding  to  the  second  and  third  transmis¬ 
sion  resonance. 


4.3  Conclusions 

We  have  observed  TH  signal  from  arrays  of  sub- wavelength  apertures.  The  peak  TH  is  obtained  at 
incidence  angles  corresponding  to  EOT  of  the  fundamental,  with  signal  strength  comparable  to  SH 
under  strong  symmetry-breaking  conditions.  There  is  a  clear  effect  of  lattice  spacing  and  aperture 
size  on  TH  signal,  where  an  anomaly  is  observed  in  the  scaling  of  TH  with  aperture  size  at  different 
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incidence  angles  corresponding  to  EOT. 

Future  work  will  focus  on  multi-photon  luminescence,  the  prospects  for  super-continuum  gener¬ 
ation,  and  the  incorporation  of  strongly  nonlinear  materials  within  the  apertures,  such  as  quantum 
dots. 


20 


References 

[1]  V.  Govindan  and  S.  Blair  “Limitations  on  nonlinear  pulse  propagation  in  coupled-resonator 
slow-light  waveguides,”  Optics  Express  15,  3922-3930  (2007). 

[2]  V.  Govindan  and  S.  Blair  “Nonlinear  pulse  interaction  in  microresonator  slow- light  waveg¬ 
uides,”  Journal  of  the  Optical  Society  of  America  B  25,  C23-C30  (2008). 

[3]  V.  Govindan  and  S.  Blair  “Analysis  of  optical  arma  architectures  in  the  slow-light  regime,” 
Journal  of  the  Optical  Society  of  America  B  25,  C116-C126  (2008). 

[4]  J.  Goeckeritz  and  S.  Blair  “One-dimensional  photonic  crystal  rib  waveguides,”  Lightwave  Tech¬ 
nology,  Journal  of  25,  2435-2439  (2007). 

[5]  J.  Goeckeritz  and  S.  Blair  “Optical  characterization  of  coupled  resonator  slow-light  rib  waveg¬ 
uides,”  Optics  Express  18,  18190-18199  (2010). 

[6]  E.  Popov,  M.  Neviere,  J.  Wenger,  P.-E.  Lenne,  H.  Rigneault,  P.  Chaumet,  N.  Bonod, 
J.  Dintinger,  and  T.  Ebbesen  “Eield  enhancement  in  single  subwavelength  apertures,”  Journal 
of  the  Optical  Society  of  America  A  23,  2342-2348  (2006). 

[7]  F.  Mahdavi,  Y.  Liu,  and  S.  Blair  “Modeling  fluorescence  enhancement  from  metallic  nanocav¬ 
ities,”  Plasmonics  2,  129-141  (2007). 

[8]  T.  Xu,  X.  Jiao,  G.  P.  Zhang,  and  S.  Blair  “Second-harmonic  emission  from  sub-wavelength 
apertures;  Effects  of  aperture  symmetry  and  lattice  arrangement,”  Optics  Express  15,  13894- 
13906  (2007). 

[9]  T.  Xu,  X.  Jiao,  and  S.  Blair  “Third-harmonic  generation  from  arrays  of  sub-wavelength  metal 
apertures,”  Opt.  Express  17,  23582-23588  (2009). 

[10]  L.  V.  Hau,  S.  E.  Harris,  Z.  Dutton,  and  C.  H.  Behroozi  “Light  speed  reduction  to  17  metres 
per  second  in  an  ultracold  atomic  gas,”  Nature  397,  594-598  (1999). 

[11]  A.  Yariv,  Y.  Xu,  R.  K.  Lee,  and  A.  Scherer  “Coupled-resonator  optical  waveguides:  a  proposal 
and  analysis,”  Optics  Letters  24,  711-713  (1999). 

[12]  J.  Heebner,  R.  W.  Boyd,  and  Q.-H.  Park  “SCISSOR  solitons  and  other  novel  propagation 
effects  in  microresonator-modihed  waveguides,”  Journal  of  the  Optical  Society  of  America  B 
19,  722-731  (2002). 

[13]  A.  Melloni,  F.  Morichetti,  and  M.  Martinelli  “Linear  and  nonlinear  pulse  propagation  in  cou¬ 
pled  resonator  slow-wave  optical  structures,”  Optical  and  Quantum  Electronics  35,  365-379 
(2003). 

[14]  Y.  Chen  and  S.  Blair  “Nonlinearity  enhancement  in  hnite  coupled-resonator  slow-light  waveg¬ 
uides,”  Optics  Express  12,  3353-3366  (2004). 

[15]  M.  A.  Allen  and  G.  S.  Kino  “On  the  theory  of  strongly  coupled  cavity  chains,”  IEEE  Trans¬ 
actions  on  Microwave  Theory  and  Techniques  8,  362-372  (1960). 

[16]  H.  J.  Curnow  “A  general  equivalent  circuit  for  coupled-cavity  slow-wave  structures,”  IEEE 
Transactions  on  Microwave  Theory  and  Techniques  13,  671-675  (1965). 


21 


[17]  L.  Maleki,  A.  B.  Matsko,  A.  A.  Savchenkov,  and  V.  S.  Ilchenko  “Tunable  delay  line  with 
interacting  whispering-gallery- mode  resonators,”  Optics  Letters  29,  626-628  (2004). 

[18]  R.  W.  Boyd,  D.  J.  Gauthier,  A.  L.  Gaeta,  and  A.  E.  Willner  “Maximum  time  delay  achievable 
on  propagation  through  a  slow-light  medium,”  Physical  Review  A  71,  (2005). 

[19]  R.  S.  Tucker,  P.-G.  Ku,  and  C.  J.  Ghang-Hasnain  “Slow-light  optical  buffers;  Gapabilities  and 
fundamental  limitations,”  Journal  of  Lightwave  Technology  23,  4046-  (2005). 

[20]  J.  B.  Khurgin  “Optical  buffers  based  on  slow  light  in  electromagnetically  induced  transparent 
media  and  coupled  resonator  structures:  comparative  analysis,”  Journal  of  the  Optical  Society 
of  America  B  22,  1062-1074  (2005). 

[21]  F.  G.  Sedgwick,  G.  J.  Ghang-Hasnain,  P.  C.  Ku,  and  R.  S.  Tucker  “Storage-bit-rate  product 
in  slow-light  optical  buffers,”  Electronics  Letters  41,  1347-1348  (2005). 

[22]  Y.  Ghen,  G.  Pasrija,  B.  Farhang-Boroujeny,  and  S.  Blair  “Engineering  the  nonlinear  phase  shift 
using  multi-stage  auto-regressive  moving- average  optical  filters,”  Applied  Optics  13,  2564-2574 
(2005). 

[23]  G.  Pasrija,  Y.  Chen,  B.  Farhang-Boroujeny,  and  S.  Blair  “DSP  approach  to  the  design  of 
nonlinear  optical  devices,”  EURASIP  Journal  on  Applied  Signal  Processing  2005,  1485-1497 
(2005). 

[24]  C.  K.  Madsen  and  J.  H.  Zhao  Optical  Filter  Design  and  Analysis:  A  Signal  Processing  Ap¬ 
proach.  Wiley  1999. 

[25]  S.  Darmawan,  Y.  Landobasa,  and  M.  Chin  “Phase  engineering  for  ring  enhanced  mach-zehnder 
interferometers,”  Optics  Express  13,  4580-4588  (2005). 

[26]  E.  M.  Dowling  and  D.  L.  MacFarlane  “Lightwave  lattice  filters  for  optically  multiplexed  com¬ 
munication  systems,”  Lightwave  Technology,  Journal  of  12,  471-486  (1994). 

[27]  K.  Jinguji  “Synthesis  of  coherent  two-port  optical  delay  line  circuits  with  ring  waveguides,” 
Journal  of  Lightwave  Technology  14,  1882-1898  (1996). 

[28]  A.  Mellon!  and  M.  Martinelli  “Synthesis  of  direct-coupled-resonators  bandpass  filters  for  WDM 
systems,”  Journal  of  Lightwave  Technology  20,  296-303  (2002). 

[29]  T.  W.  Ebbesen,  H.  J.  Lezec,  H.  F.  Ghaemi,  T.  Thio,  and  P.  A.  Wolff  “Extraordinary  optical 
transmission  through  sub-wavelength  hole  arrays,”  Nature  391,  667-669  (1998). 

[30]  L.  Salomon,  F.  Grillot,  A.  V.  Zayats,  and  F.  Fornel  “Near-held  distribution  of  optical  transmis¬ 
sion  of  periodic  subwavelength  holes  in  a  metal  him,”  Physical  Review  Letters  86,  1110-1113 
(2001). 

[31]  Y.  Liu  and  S.  Blair  “Fluorescence  enhancement  from  an  array  of  sub-wavelength  metal  aper¬ 
tures,”  Optics  Letters  28,  507-509  (2003). 

[32]  A.  Nahata,  R.  A.  Linke,  T.  Ishi,  and  K.  Ohashi  “Enhanced  nonlinear  optical  conversion  using 
periodically  nanostructured  metal  hlms,”  Optics  Letters  28,  423-425  (2003). 


22 


[33]  N.  Rakov,  F.  E.  Ramos,  and  M.  Xiao  “Strong  second-harmonic  generation  from  a  thin  silver 
film  with  randomly  distributed  small  holes,”  Journal  of  Physics:  Condensed  Matter  15,  L349- 
L352  (2003). 

[34]  M.  Airola,  Y.  Liu,  and  S.  Blair  “Second-harmonic  generation  from  an  array  of  sub-wavelength 
metal  apertures,”  Journal  of  Optics  A  7,  S118-S123  (2005). 

[35]  W.  Fan,  S.  Zhang,  N.  C.  Panoiu,  A.  Abdenour,  S.  Krishna,  J.  R.  M.  Osgood,  K.  J.  Malloy, 
and  S.  R.  J.  Brueck  “Second- harmonic  generation  from  a  nanopatterned  isotropic  nonlinear 
material,”  Nano  Letters  6,  1027-1030  (2006). 

[36]  A.  Lesuffleur,  L.  K.  S.  Kumar,  and  R.  Gordon  “Enhanced  second  harmonic  generation  from 
nanoscale  double-hole  arrays  in  a  gold  film,”  Applied  Physics  Letters  88,  261104  (2006). 

[37]  J.  A.  H.  van  Nieuwstadt,  M.  Sandtke,  R.  H.  Harmsen,  F.  B.  Segerink,  J.  C.  Prangsma,  S.  Enoch, 
and  L.  Kuipers  “Strong  modification  of  the  nonlinear  optical  response  of  metallic  subwave¬ 
length  hole  arrays,”  Physical  Review  Letters  97,  146102  (2006). 

[38]  A.  Lesuffleur,  L.  K.  S.  Kumar,  and  R.  Gordon  “Apex-enhanced  second  harmonic  generation 
by  using  double- hole  arrays  in  a  gold  film,”  Physical  Review  B  75,  045423  (2007). 

[39]  P.  Miihlschlegel,  H.-J.  Eisler,  O.  J.  F.  Martin,  B.  Hecht,  and  D.  W.  Pohl  “Resonant  optical 
antennas,”  Science  308,  1607-1609  (2005). 

[40]  S.  Kim,  J.  Jin,  Y.-J.  Kim,  I.-Y.  Park,  Y.  Kim,  and  S.-W.  Kim  “High-harmonic  generation  by 
resonant  plasmon  field  enhancement,”  Nature  453,  757-760  (2008). 

[41]  X.  Jiao,  J.  Goeckeritz,  S.  Blair,  and  M.  Oldham  “Localization  of  near-held  resonances  in 
bow-tie  antennae:  inhuence  of  adhesion  layers,”  Plasmonics  4,  37-50  (2009). 

[42]  H.  Aouani,  J.  Wenger,  D.  Gerard,  H.  Rigneault,  E.  Devaux,  T.  W.  Ebbesen,  F.  Mahdavi,  T.  Xu, 
and  S.  Blair  “Grucial  role  of  the  adhesion  layer  on  the  plasmonic  fluorescence  enhancement,” 
AGS  Nano  3,  2043-2048  (2009). 


23 


